Hydrogen-bonded Silica Gels Dispersed in a Smectic Liquid Crystal: A Random Field 

XY System 
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The effect on the nematic to smectic-A transition in octyicyanobiphenyi (8CB) due to dispersions 
of hydrogen-bonded siiica (aerosii) particies is characterized with high-resoiution x-ray scattering. 
The particies form weak gels in 8CB creating a quenched disorder that replaces the transition with 
the growth of short range smectic correlations. The correlations include thermal critical fluctuations 
that dominate at high temperatures and a second contribution that quantitatively matches the static 
fluctuations of a random fleld system and becomes important at low temperatures. 
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PACS numbers: 64.70.Md, 61.30.Eb, 61.10.-i 



Quenched disorder in condensed matter systems can 
profoundly affect their thermodynamic behavior, creat- 
ing unique challenges for statistical mechanics. For flu- 
ids, confinement within a random porous medium offers 
a convenient method for introducing such disorder, and 
studies of superfluid helium EJ and binary liquidsBI in aero- 
gels - rigid, highly-porous chemically-bonded silica gels 
- have led to insights into the implications of disorder on 
phase transitions. Similarly, the fragile nature of the me- 
somorphic phases of liquid crystals and their sensitivity 
to surface interactions make these systems well suited for 
studies within porous media. In particular, the nematic 
(N) to smectic-A (SmA) transition, which typically ex- 
hibits anisotropic 3D-XY critical behavior, offers an ideal 
case for understanding the effects of a random environ- 
ment on continuous symmetry fluctnations. The dernfln- 
stration through x-ray scatteringBu and calorimetryd of 
the destruction of the N-SmA transition by confinement 
in aerogels has shed light on these effects. Detailed the- 
oretical work has further elucidated the fragility of the 
smectic phase to arbitrarily weak quenched disorder and 
has introduced the possibility that for sufficiently weak 
disorder novel "Bragg glass" phases may become stableQ. 

In an effort to probe liquid crystals with increasingly 
weak disorder in a controlled way, we have conducted 
a high-resolution x-ray scattering study of the effect of 
hydrogen-bonded silica (aerosii) gels on the N-SmA tran- 
sition in octyicyanobiphenyi (8CB). We find that the 
constraints imposed by these weakly-connected gels pro- 
foundly alter the smectic ordering, leading to behavior 
quantitatively like that induced by random fields. Ran- 
dom field systems have been fruitful models for exploring 
theoretically effects of quenched disorder, and random 
field Ising magnets have provided the primary experi- 
mental testing ground for these ideasQ. 8CB in porous 
silica thus introduces an opportunity to make experi- 
mental contact with random field theories for a transi- 



tion that breaks a continuous symmetrj^lHil3. Consis- 
tent with theoretical predictions, we find that the smec- 
tic phase is destroyed by the disorder and is replaced by 
the growth of short-range smectic correlations. The low- 
temperature correlation length exceeds the typical void 
size in the gels by approximately a factor of four, con- 
sistent with the picture of the aerosii creating a weak to 
intermediate disorder. The temperature dependence of 
the correlation function for different densities of gel fur- 
ther indicates that the disorder suppresses nematic fluc- 
tuations in the liquid crystal. 

The smectic state is distinguished from the nematic 
by a one-dimensional modulation in density, and a two- 
component order parameter specifying the amplitude and 
phase of the modulation characterizes the N-SmA tran- 
sition. To introduce quenched disorder, we disperse in 
the 8CB hydrophilic aerosii, 70 A diameter silica par- 
ticles coated with hydroxyl groups. At volume frac- 
tions above roughly 1%, these particles form a hydrogen- 
bonded thixotropic gel, and interactions between the sur- 
faces of the gel and liquid crystal molecules ( "anchoring" ) 
introduce locally preferred orientations and phases to the 
smectic ordering that compete with the globally uniform 
order of the pure smectic. Aerosii gels have strong simi- 
larities with aerogels but represent a weaker disorder for 
the smectic in two essential respects. First, the ability 
to grow the aerosii gel directly within the liquid crystal 
permits structures of higher porosity. Also, as discussed 
below, for a given silica density the hydrogen-bonded 
aerosils perturb the smectic ordering less severely than 
rigjidj-aerogels. However, consistent with recent stud- 
iegli't3 we find that aerosils above 1% volume fraction 
create a reproducible, density-dependent quenched dis- 
order. In our experiments we have tuned the strength 
of this disorder by varying the density of silica, ps , from 
0.025 grams of Si02 per cm'^ of 8CB, just above the gela- 
tion threshold, to 0.341 grams of Si02 per cm'^ of 8CB. 
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ReferenceEi provides details of the sample preparation 
and fractal structure of aerosil gels. 

To characterize the effect of the aerosil gels on the N- 
SmA transition, we have measured with x-ray scattering 
the static structure factor, S{q), which provides direct 
information about the correlation function of the smec- 
tic order. The experiments were conducted on the X20 
beamline of the National Synchrotron Light Source using 
8 keV photons. 1 mm thick samples were held between 
kapton sheets in an oven containing dry nitrogen gas, 
and the scattering was performed in transmission geome- 
try. Following prolonged exposure to the x-rays, the 8CB 
with dispersed aerosil displayed evidence of damage, in- 
cluding progressively degraded smectic order. Therefore, 
we carefully limited the exposures to levels below which 
any damage was detectable. 

Figure 1 shows representative x-ray scattering line- 
shapes for two values of ps- In accord theoretical 
predictions of Radzihovsky, Toner, et alQx3, for all den- 
sities and for the full experimental temperature range, 
the diffraction peak corresponding to smectic layering is 
markedly broader than the resolution-limited peak mea- 
sured for pure 8CB in the SmA phase, demonstrating the 
destruction of the smectic long range order by the aerosil. 
We model this short-ranged order with a structure factor 
containing two contributions: 
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where the wave vector q|| is along the smectic layer nor- 
mal and q± is perpendicular to it. 

The first term in S{q), an anisotropic Lorentzian with 
quartic corrections, is the form for the susceptibility in 
pure nematic 8CB and describes the critical therpial fluc- 
tuations on approaching the N-SmA transitionlU. The 
second term, which is simply proportional to this suscep- 
tibility squared, is designed to account for static fluctua- 
tions induced by the quenched disorder. This term is mo- 
tivated by studies on random field Ising magnets in which 
such an expression has been shown to describe accurately 
the shoijt^range correlations induced by the st|atiGjfluc- 
tuationsta and has been justified theoreticallyll30. In 
the limit of long range order, the second term evolves 
into a Bragg-peak, <72S{q\\ — qa)d{q±_'p3. The soUd lines 
in Fig. 1 are the results of fits to the powder average 
of Eq. (1) convolved with the experimental resolution 
function. The measured scattering profile also includes 
a temperature-independent, sloping background from the 
aerosil gel, shown in the figures as dashed lines. As Fig. 1 
illustrates, Eq. (1) provides an excellent description of the 
lineshape over a wide range of temperatures and aerosil 
densities. In these fits, we have taken and c as func- 
tions of ^||, with and c(5||) set by theirrelations 
in pure 8CB, which are known to high precisionEJ. Thus, 
each fit has only four free parameters - qo, £,\\, cti, and 



a2 - with two parameters, ^|| and (72 /ai, controlling the 
shape of the peaks. 

At temperatures near and above the N-SnxA. transition 
temperature for pure 8CB, T^^ = 306.98A'y, the mea- 
sured lineshapes in the aerosil samples are described very 
well by the thermal fluctuation term alone (i.e., (72 = 0). 
Thus, for T > T^^, the smectic correlations in aerosil 
samples are very similar to the critical fluctuations of the 
pure liquid crystal. Figure 2(a) displays ^|| for a series 
of aerosil densities in this high temperature region. The 
solid line in the figure, ^|| for pure 8CB (aligned in the 
nematic by a magnetic field)Lj, demonstrates the strong 
similarity between thermal fiuctuations in the pure liquid 
crystal and in the liquid crystal with dispersed aerosil. 
At these high temperatures, ^|| is smaller than the av- 
erage void size, lo, in the gels, and the infiuence of the 
gel is minimal. (The values of lo range, from 200 A for 
Ps = 0.341 to 2650 A for ps = 0.025EI.) Rather than 
follow the diverging behavior associated with a thermo- 
dynamic transition, however, for 8CB with dispersed 
aerosil deviates from that of pure 8CB as the tempera- 
ture decreases, and the smectic correlations remain short- 
ranged. Figure 2(b) shows the behavior of for four 
values of ps to low temperatures. The susceptibility, cti, 
for 8CB with dispersed aerosil similarly tracks that of the 
pure system at high temperatures but fails to diverge. 

We stress the importance of including the quartic term 
(c > 0) in Eq. (1) to account correctly for the tempera- 
ture dependence of the lineshapes. The values of c vary 
from approximately 0.25 for ^y = 50A, its value 4 K 
above T°^, to less than 0.02 for C|| > IO^'AEI. At high 
temperatures, when c ~ 0.25, the thermal fiuctuation 
term is effectively the square of a Lorentzian along q±, 
leading to a powder-averaged lineshape that strongly re- 
sembles a Lorentzian. As increases on cooling and c 
decreases, the powder average of the first term in Eq. (1) 
deviates strongly from a Lorentzian, and the tails of the 
lineshape become more pronounced. We note that the x- 
ray lineshapes of 8CB in rigid aerogels resemble this be- 
havioila. The analysis by Rappaport et al. of the aerogel 
results, however, differs significantly from that presented 
here. In particular, no measurable signal fropji thermal 
fiuctuations was reported for 8CB in aerogelHa. 

At temperatures below T^^, deviations from the form 
describing the critical behavior of pure 8CB become ap- 
parent in the lineshape of 8CB with dispersed aerosil, 
particularly for large ps- These deviations indicate an in- 
creasing importance of interactions with the gel, rather 
than intrinsic thermal fluctuations, in determining the 
smectic correlations. The second term in Eq. (1) describ- 
ing random-fleld induced, short-range order accounts for 
these changes in the correlations very well. Figure 3 dis- 
plays f72 as a function of temperature for two represen- 
tative values of ps- The rapid growth of f72 on cooling 
reflects the dominance of these random-field static fluctu- 
ations at low temperatures. In this low temperature re- 
gion, reaches an essentially temperature-independent 
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value, as Fig. 2(b) illustrates. Figure 4 displays the low 
temperature correlation lengths, ^jj'"^ as a function of ps- 
The solid line in Fig. 4 is the result of a fit to the nematic 
correlation lengths of the homolog £CB with dispersed 
aerosils measured by Bellini et a/£j with static light 
scattering. These light scattering studies demonstrate 
that the nematic order in a liquid crystal with dispersed 
aerosil breaks up into very large, but finite, domains. 
Since the smectic layering forms within these nematic 
domains, the nematic domain sizes naturally set upper 
limits for the range of the smectic correlations. However, 
S^^'^ remains well below these limits for all ps, demon- 
strating the sensitivity of the smectic phase to quenched 
disorder. The dashed line in Eig. 4 is the mean aerosil 
void size, /q, as a function of pjL5. For the range of densi- 
ties studied, the ratio of ^jj'"^ to Iq is nearly independent 
of density and is roughly equal to four, consistent with an 
intermediate strength disordering field. Tiiis result con- 
trasts with the smectic ordering in aerogelB for which the 
measured correlation lengths barely exceed the void sizes, 
illustrating again the weaker nature of disorder induced 
by aerosil. 

At low temperatures (where c approaches a constant 
value) cr2 is proportional to the integrated intensity of 
the static fluctuation term in <5'(g). The growth of this 
intensity strongly resembles that of an order parameter 
squared, and the solid lines in Fig. 3 are the results of fits 
to the form (T2 oc (T2 — T)^, where both x and T2 depend 
on Ps . The inset of Fig. 3 displays the change in x with ps . 
We emphasize that x is an effective exponent since the 
aerosil destroys the true N-SmA phase transition. How- 
ever, associating the integrated intensity from the static, 
short-range correlations with an order parameter squared 
(so that X — 2/3), we observe that the effective exponent 
spans the range from near the tricritical value (like with 
pure 8CB) at small ps to 3D XY-like at large ps- For 
pure liquid crystals, the measured critical behavior at 
the N-SmA transition is affected by couplings between 
the nematic and smectic order parameters (de Gennes 
coupling)E3. The variation of x with ps suggests that the 
quenched disorder suppresses the nematic susceptibility 
that drives the critical behavior of the pure system away 
from 3D-XY universality. This result agrees nicely with 
specific heat studies of 8CB with dispersed aerosil which 
show the effective heat capacity exponent, o^jdecreasing 
toward the 3D-XY value with increasing psIiilO. 

In conclusion, these x-ray scattering studies have 
shown that 8CB with dispersed aerosil provides a model 
random field system for a transition that breaks a con- 
tinuous symmetry. In particular, the success of Eq. (1) in 
describing the smectic correlations and the corresponding 
behaviors of ^|| and (T2 form a compelling picture in which 
the aerosil gel creates a random local-orienting field that 
destroys the N^SmA transition, consistent with theoreti- 
cal predictionaS. In the random field Ising magnet, slow 
dynamics and,viQn-equilibrium effects strongly influence 
the behaviorc3'E2l. Efforts to observe hysteretic or time 



dependent behavior indicative of metastability in these, 
scattering measurements or in the calorimetry studiesEll 
on 8CB with dispersed aerosil revealed no such effects. 
These results suggest that the continuous symmetry al- 
lows the system to find equilibrium efficiently, even in the 
presence of random fields, in marked contrast to the be- 
havior in random field Ising systemsQ. However, a better 
understanding of the dynamics would be beneficial for 
placing the observed ordering in the context of equilib- 
rium statistical mechanics. 
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FIG. 1. X-ray scattering intensity for 8CB with dispersed 
aerosil for = 0.025 g/cm^ at (a) 288.5 K and (b) 307.6 
K and = 0.282 g/cm^ at (c) 300.6 K and (d) 306.6 K. 
The solid lines are the results of fits to the powder average of 
Eq. (1) convolved with the experimental resolution and the 
dashed lines are the temperature independent background due 
to scattering from the aerosil gel. The higher temperature 
lineshapes, (b) and (d), are due to thermal critical fluctua- 
tions as in pure 8CB while the low temperature lineshapes, 
(a) and (c), have forms expected for short-range correlations 
dominated by static, random-field induced fluctuations. 

FIG. 2. The parallel correlation length, ^y, as a function 
of temperature (a) at high temperature for a series of aerosil 
densities: p, = 0.025 (Q), 0.041 (•), 0.051 (filled diamond), 
0.078 (open box), 0.105 (V), 0.161 (A), 0.220 (filled box), 
0.282 (filled triangle), and 0.341 (filled upside-down triangle); 
and (b) on an expanded scale to low temperature for four 
values of ps- The solid line is for pure 8CB [14]. 

FIG. 3. The integrated intensity of the static fluctuation 
term, a2, as a function of temperature for ps = 0.025 (•) and 
Ps = 0.220 (O)- The absolute magnitudes of a2 for each ps 
are arbitrary and should not be compared. The solid lines 
arc fits to (T2 cx (T2 — T)^. The inset shows x versus ps. The 
effective exponent x spans the range for an order parameter 
exponent from near the tricritical value, 0.5, at small ps to 
the 3D XY value, 0.69, at large ps- 



FIG. 4. Low temperature correlation lengths, versus 
Ps- The solid line is the nematic correlation length (for 6CB) 
[19], and the dashed line is the mean aerosil void size, Iq. 
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Figure 1 S. Park et al. 




Figure 2 S. Park et al. 
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Figure 3 S. Park et al. 
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Figure 4 S. Park et al. 




